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We investigate the possibilities of generation of solitons 
and vortices in a degenerate gas of neutral fermionic atoms. 
In analogy with, already experimentally demonstrated, tech- 
nique applied to gaseous Bose-Einstein condensate we propose 
the phase engineering of a Fermi gas as a practical route to 
excited states with solitons and vortices. We stress that soli- 
tons and vortices appear even in a noninteracting fermionic 
gas. For solitons, in a system with sufficiently large number of 
fermions and appropriate trap configuration, the Pauli block- 
ing acts as the interaction between particles. 

PACS number(s): 05.30.Fk, 03.75.Fi 

Recent experimental achievement of quantum degen- 
eracy in a dilute gas of fermionic 40 K atoms [|l] has trig- 
gered theoretical interest in properties of such systems. 
In particular, the interactions between ultracold fermions 
were studied both in the static and dynamic context. 
Since Fermi atoms in the same spin state do not inter- 
act via s-wave collisions, having atoms in two hypcrfinc 
states is necessary to perform cooling of a system us- 
ing evaporative technique. That is the idea implemented 
in experiment Q. In Ref. |2j in phase collective excita- 
tions of such two-component system are discussed within 
the hydrodynamic limit and method for detection of the 
onset of Fermi degeneracy is proposed. Static proper- 
ties of two-component system were already analyzed in 
Ref. ||. In a one-component gas in the absence of s- 
wave collisions other types of forces are getting impor- 
tant, for example, the dipole-dipole forces. The proper- 
ties of trapped fermionic dipoles, including the stability 
analysis, have been investigated recently in Ref. In- 
teresting results on the transition temperature to BCS 
phase and the detection of Cooper pairing were also ob- 
tained B. 

In this Letter, we investigate another aspect of dy- 
namic behavior of a degenerate Fermi gas. Following re- 
cent success in experimental demonstration of topological 
defects in Bose-Einstein condensate such as vortices p|]7j 
we ask the question whether vortices can be observed in 
a gas of neutral fermionic atoms. We are interested in 
a generation of solitons and vortices in a normal state 
of a Fermi gas, i.e., above the temperature for a BCS- 
type phase transition. We find that the phase imprinting 
technique [^H already applied for the Bose-Einstein con- 
densate also works for an ultracold Fermi gas, although 
the nature of vortices appears to be different. In fact, 



because of technical reason only dark solitons were ob- 
served in the Bose-Einstein condensate so far by using 
this method. Then, we analyze the "birth" and dynam- 
ics of solitons in a fermionic gas too and discover some 
differences in comparison with bosons. 

As opposed to the Bose-Einstein condensate, the 
fermionic system possesses no macroscopic wave func- 
tion. This fact does not exclude, however, the existence 
of excited states of such systems with solitons and vor- 
tices. To proof that, we have solved the many-body 
Schrodinger equation for a one-dimensional and three- 
dimensional noninteracting Fermi atoms in a harmonic 
trap. At zero temperature the many-body wave func- 
tion is given by the Slater determinant with the lowest 
available one-particle orbitals occupied. By using the 
sufficiently fast phase imprinting technique each atom 
acquires the same phase (j)(r), hence the orthogonality of 
one-particle orbitals is not broken. Time (unitary) evolu- 
tion does not spoil this orthogonality too, so the diagonal 
part of one-particle density matrix is always the sum of 
one-particle orbitals densities during the evolution. 

In one-dimensional case only generation of solitons can 
be considered. We take the phase being imprinted in the 
form 4>(z) = 0o(l+tanh(z/£))/2 with the phase jump de- 
termined by 4>q and the width of the phase jump equal to 
£. In Fig. [| we plotted the density profiles of a system of 
N=400 noninteracting Fermi atoms at various times. It is 
clear that two solitons, the dark and the bright one, are 
generated. They propagate in opposite directions with 
the speed of bright soliton being higher. After reaching 
the edge of the system the solitons return to the center 
of a trap and oscillatory motion of solitons is observed 
(oscillatory motion of solitons in Bose-Einstein conden- 
sate was reported in Ref. E0]). It is important to note 
that in each graph of Fig. |l[we show two curves. One of 
them is a result of a solution of many-body Schrodinger 
equation. The second comes from a one-dimensional 
Thomas-Fermi approximation (discussed below). The 
perfect agreement between an exact and approximate so- 
lutions certainly proves the validity of the Thomas-Fermi 
approach in one dimension for sufficiently large number 
of atoms. It also proves the existence of solitons in the 
system of noninteracting fermions suggested recently in 
Ref. Jll]] . It means also that for N=400 atoms and more, 
the Pauli blocking plays the role of the interaction be- 
tween atomic fermions and assures the local equilibrium 
necessary within the Thomas-Fermi approximation. 

Since the harmonic potential is a special one in a sense 
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that it allows for a propagation of wave packets without 
spreading, we have also checked the case of a system of 
fermions confined in a one-dimensional box. Again, af- 
ter phase imprinting we found two well-defined pulses 
propagating in opposite directions with distinct veloci- 
ties. The speed of the bright soliton is higher than the 
speed of sound in a fermionic gas at T=0 temperature 
(given by (Tm/m) p) whereas the dark soliton propagates 
slower than the sound wave. 
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FIG. 1. Density profiles of ID noninteracting fermionic gas 
at different times: (a) 0.16, (b) 0.31, (c) 0.47, and (d) 0.63 in 
units of 1/tJ after imprinting a single phase step of 4>o = 2.0 it 
and £ = 0.5 osc. units. Note that each graph consists ac- 
tually of two curves: one obtained based on the many-body 
Schrodinger equation (solid line) and the second coming from 
the one-dimensional Thomas-Fermi model (dashed line). 

It is important to verify the above scenario in more 
than one-dimensional space. To this end, let's con- 
sider a noninteracting Fermi gas in an arbitrary three- 
dimensional harmonic trap characterized by the frequen- 
cies Uj, u) y , and lj z . One-particle orbitals are taken as 

a product cpn} (x) (pny (y) <fn} (z) of eigenvectors of the 
Hamiltonians of one-dimensional harmonic oscillators. 
The phase is being imprinted in the same form as in 
one-dimensional case. Since the Hamiltonian separates 
in coordinates 'x', 'y', and V, the time propagation of 
each orbital after phase imprinting can be easily reduced 
to one-dimensional problem: 

where 

and H z is the Hamiltonian of one-dimensional harmonic 



oscillator. Now the diagonal part of one-particle density 
matrix is given by the expression: 

^*) = 4 £ b&M)! 3 |<>M)| a 

E„ x » y „ 2 <E F 

xkSOM)! 2 (i) 

where N is the number of atoms and sum is performed 
over the one-particle states below Fermi energy. 

One of the obvious regimes where the Thomas-Fermi 
approximation holds is that of a cigar-shaped trap elon- 
gated enough, depending on the number of confined 
atoms. For example, in the case of a ground state of 
N=400 fermions in the harmonic potential with oscil- 
lator frequencies taken as uj x = u) y = AQQlo z only the 
smallest possible radial quantum numbers are involved: 
n x = n y = 0. It follows from the formula ([!]) that 
the density along V axis (or any line parallel to it) is 
just one-dimensional density multiplied by a constant. 
Certainly, the overall behavior found in one-dimensional 
case translates to three-dimensional one. Conditions pre- 
sented above could be considered also as a criterion for 
treating the real Fermi gas as a one-dimensional system. 
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FIG. 2. Images of the fermionic density of N=1000 atoms 
after writing a single phase step of 1.0 7r and £ = 0.2 osc. units 
onto a three-dimensional noninteracting Fermi gas confined 
in a trap configured as follows: u x = 2o; z and uj v = 64o; z . 
The snapshots are taken at times of (a) 0.10, (b) 0.25, (c) 
0.40, and (d) 0.60 in units of l/u a . 

Now we start to depart from "one-dimensional" ge- 
ometry. In Figs. H and || we show the contour plots of 
density in the 'xz' plane and the density profiles along V 
direction respectively after writing a single phase step 
onto a three-dimensional gas, although confined in a 
disk-shaped trap. The system behaves effectively as a 
two-dimensional gas; the number of atoms and the fun- 
damental frequencies are chosen in such a way that one 
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{fi y ) of the quantum numbers equals zero. In this case 
we again observe two pulses moving in opposite direc- 
tions, however they are broader in comparison with one- 
dimensional case (see Figs, [j] and ||). At the same time 
results obtained based on the Thomas-Fermi approxima- 
tion differ qualitatively, showing the existence of the sec- 
ond dark soliton. Only going to much larger number of 
atoms improves the agreement between both approaches, 
however then all structures are getting less pronounced. 
In three-dimensional spherically symmetric trap we have 
never observed double pulse behavior. Instead of that, 
the dark soliton and the density wave are generated (see 
Fig. ||) much as in the case of solitons generated in the 
Bose-Einstein condensate The dark soliton propa- 
gates with the velocity smaller than the speed of sound 
and experiences (as in one-dimensional case) the oscilla- 
tory motion. 
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FIG. 3. Evolution of the density distribution of 1000 Fermi 
atoms in three-dimensional harmonic trap ( u) x — 2u z and 
uj y — 64u> z ) after imprinting a single phase step of (j>o = 1.07T 
and £ = 0.2 osc. units. The successive frames correspond to 
moments: (a) 0.25, (b) 0.40, (c) 0.60, and (d) 0.85 in units of 

1/iOx- 

A convenient way to discuss Thomas-Fermi approxi- 
mation is to start from a set of equations for reduced 
density matrices. Since there is no interaction between 
particles, the equation of motion for the one-particle den- 
sity matrix does not involve the two-particle density ma- 
trix p 2 and is given by 



+ [Vext {n , t) - V ext (r 2 ,t)]pi (n ,r 2 ,t) , 



(2) 



where V ex t(r,t) is the external potential. In the limit 
Ti — > f 2 Eq. (g) leads to the continuity equation 



dp(r, t) 

at 



V -[p{r,t)v{r,t)} =0 



(3) 



where the density and velocity fields are defined as fol- 
lows: 

p(f, t) = Urn pi (ri , r 2 , t) 



v(r,t) = — lim (Vi 
zm fx— »r 2 



V 2 )xi(ri,r2,f) 



(4) 



and xi(n-)^2,i) is the phase of the one-particle density 
matrix. 

One can also rewrite Eq. (||) introducing the center- 
of-mass (R = (n + r 2 )/2) and the relative position (s = 
fi — r 2 ) coordinates. By taking the derivative of Eq. 
(Q) with respect to the coordinate s the hydrodynamic 
Euler-type equation of motion is obtained in the limit 
s^0 



dv(r, t) 
dt 



ttt-tVh (r, t) . (5) 

mp(r,t) 2 y J m ' 



Here, the kinetic-energy stress tensor T, whose elements 
are given by 

Tki = hm — , 

m s— >o oskOsi 

depends on the one-particle density matrix (more pre- 
cisely on its amplitude oi{ri,r 2 ,t)) and is calculated 
based on a local equilibrium assumption which is the sub- 
stance of Thomas-Fermi approximation. 



o 

CO 

o 

> 

H 

CO 

z 

LU 
Q 




0.0014 
0.0012 
0.0010 
0.0008 
0.0006 
0.0004 
0.0002 



Z (osc. units) 

FIG. 4. Evolution of the density distribution of 10 4 Fermi 
atoms in three-dimensional spherically symmetric harmonic 
trap after writing a single phase step of 2.0 n and £ = 0.5 
osc. units. The successive frames correspond to moments: 
(a) 0.1, (b) 0.2, (c) 0.3, and (d) 0.4 in units of l/u. 

We have solved the set of Eqs. (||) and (^) with kinetic- 
energy stress tensor approximated by its Thomas-Fermi 
form, i.e., diagonal tensor with (in three-dimensional 
case) (?i 2 /m)(l/30/7r 2 )(67r 2 p) 5 / 3 term at each position. 
In fact, we applied the inverse Madelung transformation 
H to Eqs. © and (§) and used the Split-Operator 
method combined with the imaginary-time propagation 
technique (for generation of the ground state of the 
trapped Fermi gas) as well as the real-time propagation 
(for phase imprinting and analysis of dynamic instabili- 
ties after that). 

By using the definition (^) of the velocity field one can 
easily calculate 



3 



v(r, t) 



Yh=i Pi(r,t)vj(r,t) 



Although the one-particle velocities are irrotational, 
the global velocity is, in principle, rotational. Therefore, 
for an ultracold Fermi gas in a normal phase one should 
not expect the existence of quantized vortices. However, 
as we'll see, the excitations of fermionic gas with a hole 
in a density along some line (or a dot in two-dimensional 
case) and non-zero circulation around this line are still 
possible. 

To discuss the possibility of generation of vortices in 
a degenerate Fermi gas we involve the polar coordi- 
nates. Hence, the eigenvectors are given by: ^ nm {r, <f) — 
Ar \m\ ilr'trf) e r2 / 2 e"™^ and A 2 = n\/[n{n + \m\)l}. As 
in Refs. 0,^5[, we pass a laser pulse through the appro- 
priately tailored absorption plate before impinging it on 
the atomic gas. The laser pulse is short with a duration 
of the order of the fraction of microsecond and the fre- 
quency is detuned far from the atomic transition (only 
the dipole force is important). Under such conditions the 
atomic motion is negligible during the pulse and the only 
effect of light on the atoms is imprinting the phase. To 
generate a vortex one has to prepare the absorption plate 
with azimuthally varying absorption coefficient. Chang- 
ing the laser intensity and the duration of a pulse (the 
laser pulse area) one is able to design various vortex exci- 
tations of the Fermi gas. The state right after phase im- 
printing turns out to be dynamically unstable and soon 
after a vortex at the center of the trap is formed. In 
Fig. H we show the density after the pulse with the laser 
pulse area equal to 3 x 2ir was applied to the system of 
N=6 fermions. The stable, not quantized vortex is gener- 
ated. Increasing the number of atoms leads to the lower 
contrast of the vortex but increasing (keeping the same 
number of atoms) the pulse area restores the contrast. 
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In conclusion, we have shown that the phase imprinting 
method can be employed for generation of solitons and 
vortices in ultracold fermionic gases. We have demon- 
strated that idea explicitly by solving the many-body 
Schrodinger equation for a gas of noninteracting neutral 
fermionic atoms. There exists differences in comparison 
with the Bose-Einstein condensate case. For example, 
the presence of bright solitons moving with the veloc- 
ity higher than the speed of sound or the nature of the 
vortices which are not quantized. All these phenomena 
present for fermionic gases have nothing to do with the 
potential formation of bound states (Cooper pairs) fol- 
lowed by the condensation, since no interaction between 
Fermi atoms was included in the calculations. 
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FIG. 5. Density (in arbitrary units) of the system of N= 
fermions in the excited state with a vortex. 
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